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We observe the power law scaling behavior of the monopole 
mass and condensate in the pure compact U(l) gauge the- 
ory with the Villain action. In the Coulomb phase the 
monopole mass scales with the exponent v m — 0.49(4). In 
the confinement phase the behavior of the monopole conden- 
sate is described with remarkable accuracy by the exponent 
A:x P = 0.197(3). Possible implications of these phenomena for 
a construction of a strongly coupled continuum U(l) gauge 
theory are discussed. 



I. INTRODUCTION 

The phase transition between the confinement and 
Coulomb phases of the strongly coupled pure compact 
U(l) lattice gauge theory (compact QED) has recently 
received renewed interest and two of its aspects were 
investigated in large numerical simulations. First, sev- 
eral attempts have been made to distinguish between the 
weak first order and second order scenarios for the Wil- 
son action and in the extended coupling parameter space. 
The question is whether the two-state signal decreasing 
slowly with increasing lattice volume extrapolates to a 
nonzero or zero value, respectively, in the thermodynamic 
limit. (For a recent discussion of this subject and earlier 
references see Ref. [Q). 

Second, a scaling behavior of various bulk quantities 
and of the gauge-ball spectrum consistent with a second 
order phase transition and universality has been observed 
in the vicinity of some points on the manifold separat- 
ing the confinement and Coulomb phases, outside their 
narrow neighbourhood in which the two-state signal oc- 
curs This suggests that there may exist regions of 
the parameter space where the transition is continuous, 
though it may be not such for the particular action used 
in the simulation. 

A continuous transition would allow the construction 
of a continuum theory. But even if no critical point ex- 
ists, the theory might be considered as an effective theory, 
with finite but large cutoff, provided the range of scales 
at which a second order-like behavior holds is large. The 
question is then whether such a (possibly effective) con- 
tinuum U(l) theory would be interesting in some sense, 
e.g. would it have a phase transition, confinement, etc., 
in analogy to the lattice regularized theory. 

In this paper we address this second aspect and extend 
the investigation of the scaling behavior to observables 



related to the magnetic monopoles. To our knowledge 
this subject has not yet been investigated in a systematic 
way. But the issue is important, as whatever is interest- 
ing in the compact lattice QED is essentially related to 
the monopoles: The phase transition itself is known to 
be associated with the occurence of magnetic monopoles 
being topological excitations of the theory . Modifi- 
cations of the monopole contribution to the action have 
appreciable consequences for its position |J and prop- 
erties fin]. The long distance force in the confinement 
phase Jil| , p^| and the chiral symmetry breaking (l3| are 
best understood in terms of the monopole condensate. 
The charge renormalization in the Coulomb phase is due 
to the antiscreening by monopoles |l4|~|l6|] . Thus the exis- 
tence of an interesting effective U(l) theory presumably 
depends on whether the monopoles persist to play an 
important role in it, i.e. on the scaling behavior of the 
monopoles. 

Our findings are as follows: In the Coulomb phase 
we find at various values of the coupling (3 a very clean 
exponential decay of the monopole correlation function 
in a large range of distances. This demonstrates the 
dominance of a single particle state in this correlation 
function, the monopole, whose mass we determine. The 
monopole mass extrapolated to the infinite volume, TOoo, 
scales with the distance from the phase transition as 



where 



maaiP) = a m {/3 - /3 C ) V " 



v m = 0.49(4). 



(1.1) 



(1.2) 



This value of the exponent v m can be compared with 
the values for the correlation length exponents v obtained 
for other observables. One of these values is the non- 
Gaussian value 



0.35 



(1.3) 



found for the Lee- Yang zeros || at the transition, several 
gauge balls and approximately also for the string 
tension MM. The other, Gaussian value is 



0.5. 



(1.4) 



It has been observed previously for the scalar gauge ball 
Mm in the confinement phase. 
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Our main result is that the value (|l.2|) of v m is signifi- 



cantly greater than (1.3) and consistent with (1.4). This 
implies that monopoles would stay important in any of 
the conceivable scenarios for a construction of an effec- 
tive continuum U(l) theory. If such a theory were to be 
constructed in such a way that masses and other dimen- 
sionful observables scaling with the non-Gaussian expo- 
nent u ng were kept finite nonzero in physical units, then 
the monopoles in the Coulomb phase would get massless. 
Even if instead the Gaussian exponent v g were used, the 
monopole mass can be fixed at a finite value in physi- 
cal units. This second possibility might be particularly 
suitable in the Coulomb phase, where no other scales are 
known. 

In the confinement phase we have determined the scal- 
ing behavior of the monopole condensate extrapolated to 
the infinite volume, poc, to be 



where 



Poo = a p ((3 c - PY 



/3 cxp = 0.197(3). 



(1.5) 



(1.6) 



The function ([Lq) describes extremely well the data in a 
broad interval and the scaling behavior of the condensate 
is thus well established. 

However, the value of the magnetic exponent f3 exp 
alone is not sufficient for considering the continuum limit. 
For this purpose a renormalised condensate is needed. 
A natural procedure (like e.g. in the broken <fi 4 theory) 
would be to find a pole in the monopole correlation func- 
tion in the confinement phase. We find a contribution 
suggesting such a pole, but the data are consistent with 
its amplitude, as a function of the lattice volume, ex- 
trapolating to zero i n th e ther mod ynamic limit. Thus 
currently the results (1.5) and (L6) do not allow a con- 
clusion about the condensate in a would-be continuum 
limit. They constitute only a necessary step in this di- 
rection. 

The results are presented as follows: In Sec. 2 we sum- 
marize known facts about the Z gauge theory, which we 
use in the simulation. It is a dual equivalent to the U(l) 
lattice gauge theory with the Villain action which we ac- 
tually investigate. In the Z gauge theory the monopole 
correlation functions have a form originally found by 
Frohlich and Marchetti JlTj , which is convenient for mea- 
surements. Antiperiodic boundary conditions fl8| , |l6|| are 
used allowing the consideration of a single monopole in 
a finite volume in agreement with the Gauss law. These 
b.c. reduce the magnetic U(l) symmetry 
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to Z, 



It is pointed out that this symmetry is broken in the con- 
finement phase, which necessitates some caution during 
simulations. 

In Sec. 3 we present our results for the monopole mass 
in the Coulomb phase and determine its scaling behavior. 
In Sec. 4 the necessary extrapolation procedure of the 



monopole mass results to the infinite volume limit is dis- 
cussed. The finite size effects are sizeable, as monopoles 
with their Coulomb field are extended objects. The lead- 
ing term in the volume dependence can be determined 
from electrostatic considerations, however, and further 
terms obey a simple ansatz. 

In Sec. 5 we present the results for the monopole con- 
densate in the confinement phase and describe the search 
for a monopole condensate excitation. We discuss our re- 
sults and conclude in Sec. 6. 



II. THE DUAL FORMULATION OF PURE U(l) 
LATTICE GAUGE THEORY 

A. Z gauge theory 

The partition function of pure U(l) lattice gauge the- 
ory 



J J —IT ^ D 



(2.1) 



with an action s(6p), 9p denoting the plaquette angles, 
is related to the Z (integer) gauge theory by an exact 
duality transformation . In the Z gauge theory, the 
link variables n I(I have integer values, and the partition 
function reads 



Z= j D n exp ( — y ] s* (np) 
J P 

oo 

*>»=n e . 



xfi n XM — — oo 



with 



np = n xfJ , + n( x+ n) v - n^+c)^ — n xv 



(2.2) 
(2.3) 

(2.4) 



being the plaquette integer number associated with a pla- 
quette at position x and orientation (fJ,,v). In ( |2.2] ), 
s*(np) denotes the dual action. The dual theory is a 
gauge theory invariant under the transformations 



(V^) a 



L 



(2.5) 



with an integer valued function t x of the lattice points 
x. The theories (2.1) and (2.2) are strictly equivalent in 
the infinite volume limit and they should be comparable 
in large volumina. 

The dual action associated with the Villain action 



s{6 P ) = - log ^2 ex P 

k— — oo 



(6p + 2nky 



*t \ „2 

s {np) = 7^ n p 



(2.6) 



(2.7) 
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whereas the extended Wilson action 



C. The monopole observable 



s(9 P ) = (3 cos((9 P ) + 7 cos(2(9p) (2.8) 
corresponds to the dual action 

s*(n P ) = - log( f dz cos(z n P )e' 3 cos «+^ cos ( 2 *)) . 

J — 7T 

(2.9) 

Obviously, for a numerical simulation of the Z gauge the- 
ory the dual Villain action (|]?]) is much more practical 
than (|2.9|) . This is the reason we choose the Villain action 



(2.6) in our current work. 



B. Dual correlation functions 

The magnetic monopoles in the U(l) lattice gauge the- 
ory are described by fields <i>(x) whose correlation func- 
tions are defined by certain modifications of the partition 
function Q @ 



= \\™ exp(-^ S (#p+V P ) 



(2.10) 



Here the magnetic (in the U(l) language) flux Xp is gen- 
erated by the magnetic current density 



J 



X 



■x, 



{x+\)iiv -"-Xfiv -r -"~[x+p)vX 

~X XV \ + X^ x+ ^x^ — Xxx^ 



(2.11) 



with sources 2-7T at z\, . . . , zg and —2n at y%, . . . ,yj. For 
given sources, the existence of Xp on finite lattices de- 
pends on the topology of the lattice, i.e. on the boundary 
conditions (b.c). 



In the dual formulation, the expression (2.1C) becomes 
a usual expectation value 



= / Vn cxp^-ij^ J x ^n xl }j expf~y^s*(n.p) 



of a non-local observable 111 



exp(-»y^ J XIJ ,n x , 



(2.12) 



(2.13) 



in which J Xfl is given by the current density (2.11) 



(2.14) 



The observable in (2.12) resembles the non-local expres- 
sion of a charge operator in terms if its Coulomb field. 
In the pure U(l) gauge theory an analogous expression 
would not be gauge invariant. However, in the case of 
Z gauge theory, (2.13) is gauge invariant because the 
sources of J xtl are integer multiples of 2tt |l7| . 



To determine the correlation functions ( 2.12 ) in a 
Monte-Carlo simulation, it is useful to choose the field 
Xp so that J xfl vanishes on all links in a certain direc- 
tion t (refer red to as the "time" direction) . Then the 
observable (2.13) decomposes into factors that are local 
in t and nonlocal in the "space" directions x. 

Therefore one defines J x ^ for each source located at 
(a?o,*o) a s 

J (S,t) = 0, J(x,t) = B(x)S tt0 (2.15) 
where B{x) is a solution of 

divB(f) = 2k 8 33o (2.16) 



in three dimensions. Such a J xfi in ( [2.15 ) yields the 
monopole field 



$(£ ,*o) = exp(-i ^ J- 



Xfl'^Xfl 



(2.17) 



XfL 



Products of these fields with sources at different positions 
therefore lead to superpositions of currents J xfl in ( 2.1 2| ) 
and d2.13|). 



In the case of a source (2.15) and (2.16) the boundary 
conditions in the three space dimensions have to be cho- 



sen appropriately. Whereas there is no solution of (2.16) 
with periodic boundaries, a solution exists for antiperi- 
odic boundaries IflM 



, (x + Lej , t) = -rty (x, t) 



(2.18) 



on the lattice of spacial extension L. That solution can 
be obtained by th e thr ee-dimensional discrete Fourier- 
transformation of ( [2.16 ) on the lattice with respect to 
L-antiperio dic ba se functions. In terms of the mon opole 
observable ( 2.17 ), the boundary conditions fl2.18 ) cor- 
respond to magnetic charge conjugation C []l 6f| and are 
therefore called C-periodic. 

In principle it would be possible also for periodic b.c. 
to correlate monopole- antimonopole pairs at large rela- 
tive distance and determine approximately twice their 
mass. However, antiperiodic b.c. are much more practi- 
cal allowing to determine the mass of a single monopole 
on lattices of moderate sizes. 



D. Dispersion relation and finite size effects 

The real and imaginary part of the monopole are even 
resp. odd with respect to C. The observables with defi- 
nite momentum arc C-even 



d>+ (p, t ) = H(*(*o, to))e^° 



(2.19) 
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and C-odd 

$_(p,t )=X) 9 ( $ (^,*o))f 



(2.20) 



Since Fourier-transformation in ( 2. IE ) is done with 
L-periodic functions, the momenta are p = (p i,P2iP3) 
with pj = 2-Kkj/L and integer kj. In ( 2.20 ) Fourier- 
transformation is applied to L-antiperiodic functions, so 
that the momenta are pj = (2kj + l)n/L with integer 
kj. As a consequence, the whole monopole observable 
<f> = <f> + + <f>_ has no well-defined dispersion relation on 
finite volume lattices. Only in the infinite volume limit, 
the momentum spectra degenerate and the monopole 
mass determined by means the even and odd observables 
( 2.19 ), ( 2.20 ) lead to the same mass. We will have to take 
this fact into consideration when we examine the volume 
dependence of our results in Sec. |rvj . 

The monopole masses on finite lattices are determined 
using 



= Et 



(2.21) 



in the correlation function of the operator (2.19). In the 
case of (2.20), we assume that the mass is obtained from 
the energy by the usual particle dispersion relation 



E\ =m 2 _ +3 (2 sin— ) 



(2.22) 



However, for a particle with a Coulomb field in a finite 
volume this relation is presumably only an approxima- 
tion. We take the possible corrections into account in 
a phenomenological extrapolation to the infinite volume. 
Such a way of extrapolation to the infinite volume is nec- 
essary anyhow, as the finite size effects for a particle with 
extended Coulomb field are large and only partly under 
analytic control. The results do not change significantly if 
we use the lattice dispersion relation with 2(cosh(i?) — 1) 
instead of E 2 . 

The magnetic U(l) symmetry is reduced to a Z2 
symmetry by imposing C-periodic boundary conditions 
JUlpJ- In the confinement phase of pure U(l) gauge the- 
ory, this remaining Z2 symmetry is spontaneously bro- 
ken. This Z2 symmetry corresponds to the sign of the 
($+) expectation value whereas ) = due to C- 
antisymmetry. On finite lattices flips between both signs 
of ( < I > +} can occur, distorting the measurements. There- 
fore we discard parts of the runs where such flips occur. 



III. RESULTS IN THE COULOMB-PHASE 



A. Simulations and statistics 



We simulate the Z gauge theory with the action (2.7). 
The boundary conditions are C-periodic in spacial direc- 
tions and periodic in time direction. We use the same 



heat bath update as the authors of |16| together with a 
new implementation of the monopole observable (2.17) 
optimized for vectorization on a Cray T90. 

The lattice volumes are L 3 ■ T with T = 28 fixed and 
different L ranging from 4 to 18. Both monopole opera- 
tors ( 2.1 9| ) and (2.20) are measured after each 25 update 
sweeps for each time slice. From these data, the correla- 
tion functions are computed. In addition, we determine 
the action density. 



p ^raax Tint 



0.645 

0.647 

0.65 

0.654 

0.66 

0.668 

0.678 



18 
18 
18 
18 
16 
16 
16 



11.0 
2.5 
2.3 
1.9 
1.3 
1.1 
1.0 



TABLE I. Values of /3 at which the simulations in 
the Coulomb phase were performed. The lattice sizes are 
L = 4, 6, . . . , I/max- At each f3, at least 10 4 measurements 
have been made, Tint is the integrated autocorrelation time 
of the action density in multiples of 25 sweeps determined on 
the L = 8 lattice. 

Table | gives the simulation points in the Coulomb 
phase, the range L = 4, 6, . . . , L max of lattice sizes, and 
the integrated autocorrelation times of the action on the 
L = 8 lattice. To thermalize the system, we skipped 
the first 2.5 • 10 3 (far away from the phase transition) to 
2.5 • 10 4 sweeps (close to the phase transition). At each 
value of the parameters, at least 10 4 measurements have 
been performed. 

The restriction of lattice sizes to L < 18 is mainly due 
to the costs ~ L 6 in the determination of the monopole 
observable ( p. 17 ) because of the two summations over 
the space volume, one in the exponent, and one in the 
Fourier-transformation of ( p. 19 2.2C). The statistical er- 
rors of the monopole masses arc determined by the jack- 
knife method using 16 blocks. 



B. Determination of the monopole mass 

We use the observables and <!>_ with the low- 
est possible momenta, i.e. p = for $ + and p = 
(tt/X, n/L, tt/L) for The values of the monopole 

mass are obtained from their correlation function which 
is assumed to have approximately the form 



C±(t) = (0|*|(t)$±(0)|0> 
+ |(1|$ ± (0)|0)| 2 ( 



-ETt 



$±(0)|0)| 



(3.1) 



The first term, the monopole condensate 
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FIG. 1. Logarithmic plot of the correlation function C+ 
on the L — 12 lattice at /3 = 0.668. The solid curve is a fit 
to (O) with m+ = 0.921(5) and p = 0.0003(15). 



FIG. 2. Effective energies corresponding to the correla- 
tion function from fi gure oj. The horizontal lines show the fit 
results for from (3.1) with p = 0. 



P : 



$+|0) 



(3.2) 



is expected to vanish in the Coulomb phase in the infinite 
volume limit, but should be allowed on finite lattices. 

Figure [l] shows as an example the data obtained on 
the L = 12 lattice at j3 = 0.668 with the fit by means 



of the correlation function (3.1). The data in the whole 
range of t are consistent with a contribution of only one 
state to the correlation function. To verify this more ac- 
curately, we have determined the effective energies E e R 
with p both free and set equal to zero. The effective 
energies corresponding to the correlation function from 
figure [l] are plotted in figure |^. They are stable with 
respect to t. Similar results are obtained for both oper- 
ators ( p,19| ) and ( 2.1 9| ) at all investigated points in the 
Coulomb phase for all lattice sizes we have used. 

Thus we find that within our numerical accu racy the 
^-dependence of the correlation functions (3d) can be 
well described by a one particle contribution. This is re- 
markable, as the Coulomb phase contains massless pho- 
ton which could in principle substantially complicate the 
form of the correlation function (infra particle). It jus- 
tifies the interpretation of the observed energy as the 
energy or mass of the magnetic monopole. Furthermore, 
there is no significant deviation from p = and therefore 
from now on we quote only results obtained under the 
assumption of vanishing condensate. 



C. Scaling behaviour of the monopole mass 

From the finite volume results for the masses of the 
$ + and $_ observables, we estimate the infinite volume 
mass m-oo as described below in Sec. IV. The resulting 
values of m,oo are given in Tab. and Fig. §. 



i 

0.645 

0.647 

0.65 

0.654 

0.66 

0.668 

0.678 



0.32(4) 
0.46(2) 
0.59(3) 
0.72(2) 
0.88(3) 
1.05(3) 
1.26(2) 



TABLE II. Extrapolation m^o of the monopole mass to 
the infinite volume with ([Oj). The errors are estimated as 
described in the text. 



Figure [5] demonstrates that the monopole mass m,oc 
decreases with decreasing /3 and possibly vanishes at the 
phase transition. The smallest value of moo which we 
obtained is about moo — 1/3 (largest value of the corre- 
sponding correlation length is £ m0 n — 3). In the range 
moo = 0.3 — 1.2 a sc alin g behaviour is observed, descr ibed 
by the power law (1.1) with the critical exponent (1.2) 
and the critical point located at 



Coul 



0.6424(9). 



(3.3) 



The value oul is consistent with the result for (3 C ob- 
tained for the Villain action in HJjl|. The upper script 



"Coul" in (3.3) indicates that the position of the critical 
point has been obtained by means of an extrapolation of 
some observable from the Coulomb phase. The value of 
the critical exponent v ~ 0.5 has been found earlier to 
describe the scaling behavior of the mass of a scalar gauge 
ball in the pure U(l) lattice gauge theory with extended 
Wilson action @||. 

It is important to realize why the smallest value of moo 
we obtained is restricted. As explained in the next sec- 
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FIG. 3. Scaling of t he extrapolated monopole mass /Tic 
with the power law (1.1). 



tion, an extrapolation of the finite volume results for the 
monopole mass to the thermodynamic limit gets grad- 
ually more and more difficult as the phase transition is 
approached. In fact, if the finite volume dependence of 
the monopole masses were better understood, the phase 
transition could be further approached without entering 
the region where metastability occurs. On the 28 4 lattice 
metastability occurs only at /? < 0.6439. An extrapola- 
tion of our data using ( [L.lj) to this j3 suggests that this 
lattice size would allow to reach monopole masses at least 
as small as 0.2. Thus the range of monopole mass values 
investigated in this paper is restricted from below by an 
insufficient understanding of finite size effects and not by 
the occurence of the two-state signal. 



IV. VOLUME DEPENDENCE OF THE 
MONOPOLE MASS 

Since in the Coulomb phase there exists a long range 
interaction mediated by the massless photons, a strong 
dependence of the monopole masses on the finite spacial 
size L of the lattice is expected. Figure ||| displays the 
dependence of the masses m± obtained from the corre- 
lation functions C+(t) (solid) and C-(t) (dashed) on the 
lattices of size L. 

The mass m + which is obtained from the symmetric 
combination of monopole and antimonopole is smaller 
than 77i _. The splitting of masses apparently vanishes 
if the lattice size L gets large, as expected from the de- 
generation of the momenta of <£>+ and in the infinite 
volume. However, for smaller L the functions m±(L) 
are rather complicated, m_(L) being even nonmono- 
tonic, and currently we do not have a sufficient theoreti- 
cal understanding of this L-dependence. To extrapolate 
to L = oo we therefore combine the expected asymp- 



FIG. 4. Dependence of the monopole masses m± on 
the finite spacial lattice size L at (5 = 0.65. The solid and 
dashed cur ves for tti+ and m_, respectively, correspond to 
the ansatz (4.8) up to the l/L term. 



totic behaviour of m± (L) with various phenomenological 
ansatze. 

It seems plausible, that at larger L, where both masses 
have already similar values, the long-range Coulomb field 
of the monopole might dominate the L-dependence. So 
we try to describe m±(L) at large L by the classical en- 
ergy of a (magnetically) charged particle in a finite vol- 
ume. 

The classical energy of a (magnetic) charge g with a 
Coulomb potential (j){f) = g/Airr in a spherical volume 
of radius L in continuum is 



W cont (L) 



1 

ro 



W g = £. (4.1) 



The diverging classical energy of the point particle has 
been regularized by a restriction to r > tq ■ The classical 
consideration gives a characteristic 1 / L dependence 



m±(L) = moo - -j-- 



(4.2) 



with some value of the coefficient C\ depending on the 
shape of the finite volume, boundary conditions and reg- 
ularization. There is no difference between m + and m_ 
at the classical level. The splitting of masses is prob- 
ably a quantum mechanical effect, which we have not 
estimated. 

In order to obtain the value of ci for the cubic lattices 
used in the simulations, we first carry out a classical com- 
putation analogous to (4.1) in the lattice regularization. 
We use the lattice Coulomb potential 



^latt(^) 



exp(ip • r) 



L3 ^o 2 EU(l-cos^)' 
2tt 

Pj = —rij, rij e {0, . . . , L - 1} 



(4.3) 



on a L 3 cubic lattice with periodic boundary conditions. 
In the computation of the energy, the gradient is replaced 
by a finite difference, and the integral over the volume 
is replaced by a finite sum. We determine the energy 
numerically. The result is very well approximated by 



Wutt(L) = 0.2524-1.3881- 



2L 



(4.4) 



Compared with ( |4.1| ), only the pre- factors have been 
changed. We do not specify any errors of the coefficients 
because the statistical errors from the simulations are 
larger by an order of magnitute. The resulting estimate 
of the L-dependenc e of the monopole mass on a finite L 3 
lattice is therefore ( |4.2|) with 



ci = 1.3881 



(4.5) 



P 



Cl 



0.645 
0.647 
0.650 
0.654 
0.660 
0.668 
0.678 



0.1836 
0.1766 
0.1687 
0.1606 
0.1508 
0.1402 
0.1293 



1.3614 
1.4159 
1.4817 
1.5571 
1.6580 
1.7829 
1.9329 



1.1896 
1.2372 
1.2947 
1.3606 
1.4487 
1.5579 
1.6889 



TA BLE III. Dual renormalized coupling a mag obt aine d 
from ( |4.7| ). c\ is the coefficient (l.£) used in the ansatz (4.9). 



The drawback of this estimate is the fact that the an- 
tiperiodic boundary conditions on the gauge field are not 
fully respected. 

Another way to estimate the value of ci is to interpret 
the antiperiodic boundary conditions for the gauge field 
as a three-dimensional infinite cubic lattice of alternat- 
ing charges, the distance between nearest neighbours be- 
ing L. For the continuum space this consideration leads 
to the problem of lattice sums, relevant for various crys- 
taline materials. Taking over the well known results from 
the condensed matter physics plj implies that the energy 
of one monopol e in the field of all others is equal to the 
second term in (4.2) with 



ci = 1.7476 



*mag 



(4.6) 




P, = 0.65 



4 6 8 10 12 14 16 IB 

L 

FIG. 5. Dependence of the masses of m+ (solid line) and 
m_ (dashed li ne) o n the lattice size L at /3 = 0.65. The curves 
correspont to (4J)) up to the term 1/L 3 . 



The numerical factor in (4.6) is the Madelung constant. 
It determines the classical energy of cubic ion crystals, 
though its actual calculation requires particular mathe- 
matical attention pit]. We prefer (4.6) to (4.5) because 



it respects the correct boundary conditions for the gauge 
field. 

An estimate of a mag can be obtained from the numer- 
ical analysis of the static potential of U(l) lattice gauge 
theory with Villain action by using the Dirac relation 



g ■ e = 2n 



1 



*mag 



4a n 



(4.7) 



4| . Table |III| contains the values 
mag for different f3. 



We use the values of the renormalized electrical cou- 
pling a c i = e 2 /47r from 
of cv e i and a 

The fact that classical energy considerations suggest 
a 1/L dependence motivates our ansatze in terms of a 
polynomial in 1/L. The necessary scale is assumed to be 



provided by the infinite volume mass m c 
ansatz is 



itself. Our first 



c (±) 

m ± (L) = moo + -^j— 



„(±) 



,L 2 



m lo L3 



(4.8) 



To test the estimate (4.4), different coefficients c| + ' and 
c\ ^ are allowed even for i = 1 . Because of the theoretical 
uncertainty in this and the other ansatze, we have to be 
very cautious in estimating the errors of moo- 

Fig. H shows the fits ( Q| ) upto 1/L 2 . The results ob- 



tained for moo become stable if the fit is restricted to 
data with L > 8. In this case, moo is not signifficantly 
changed if more points at small L are omitted or if 1/ L 3 
contributions are included. The coefficients and c\ ^ 
roughly agree, but differ from the classical value c\ ob- 
tained from (4.6) and given in table III by a factor of 
1.2-1.6. 



Nevertheless, it is possible to fix c\ 
the classical value and try another ansatz 



(+) _ „(-) _ 



m±(L) 



L 



„(±) 



,L 3 



c\ to 



(4.9) 



We have used this ansatz with both estimates of the coef- 
ficient ci, given in ( p~5| ) and (4.6). It turned out that the 
results for the monopole mass in the infinite volume are 
consistent within the error bars. In the following we thus 
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If 



B = 0.645 



more reliable extrapolation procedure would be needed. 



V. RESULTS IN THE CONFINEMENT PHASE 



A. Simulations and statistics 




L 

FIG. 6. Dependence of the masses of m+ (solid line) and 
m_ (dashed line) on the lattice size L at f3 = 0.645 which is 
the point closest to the phase transition. The curves corre- 

-2 



spond to fligp up to 1/L 2 
been excluded from the fits. 



Data points with L < 10 have 



describe onl y th e results obtained using the Madclung 
constant in ( |4.6| ). 

Fig. p| s hows the same data as fig. |], but with the func- 
tions (4.9) up to 1/L 3 . The fit results become stable 
if the fit includes the 1/L 3 contributions. A restriction 
to large L is necessary only for j3 close to the phase tran- 
sition. 

We consider the values of obtained in the fit by 



means of (4.9) with all listed terms as the best determina- 
tion of the monopolc mass in the infinite volume and list 
the results in table [n]. To obtain an estimate of the errors 
of rrioo obtained in this way, we compare the fit results 
from (4.9) with and without higher 1/L k terms. Further 



we omit various number of points at small L from the 
fits. The error of moo is then estimated from the varia- 
tion of the fit results under the different conditions. The 
values of these errors are a lso given in table 0. Using 
instead of ( |4.9| ) the ansatz (4.S) results in the values of 
TOoo compatible with those in table |n| within the listed 
errors. 

Though insufficiently motivated, our extrapolation 
procedure is rather stable at the j3 values farther from 
the phas e tr ansition. E ven the use of the coefficient C\ 
given in ( 4.5 ) instead of (4.6), does not change the results 
significantly. The extrapolation is most problematic at 
the point closest to the phase transition, (3 — 0.645. Fig- 
ure |^ shows the L-dependence of the m± masses in this 
case, using (4.6). Since the finite volume mass splitting 
between m + and m_ remains significantly nonzero on 
the largest lattice, the systematic error is probably quite 
large. From the different fits we estimate it to be about 
20% at this (3. In order to reduce this error and to fur- 
ther approach the phase transition, larger lattices and a 



Our simulations in the confinement phase are per- 
formed with the same algorithm and the same boundary 
conditions as described in Sec. Ill A. The lattice volumes 
are L 3 • T with T = 24 fixed and several L. At each (3, at 
least 5000 measurements have been made starting with 
an ordered system and 5000 starting with a completely 
disordered system. Before the measurements we have 
used 2.5 • 10 3 (away from the phase transition) to 2.5 • 10 4 
(close to the transition) sweeps for thermalization. The 
monopole operators ( 2.19| ) and (2.20) are measured after 
each 25 update sweeps. 





i'max 


Tint 


0.6 


12 


1.1 


0.62 


14 


1.2 


0.625 


14 


1.3 


0.635 


14 


1.3 


0.638 


12 


2.6 


0.64 


16 


2.9 


0.641 


16 


4.2 


0.642 


16 


5.9 


0.6425 


18 


7.2 


0.643 


14 


10 


0.6435 


18 


15 


0.6436 


18 


30 


0.6437 


18 


35 



TABLE IV. Values of (3 at which the simulations in the 
confinement phase were performed. The lattice sizes are 
L = 4, 6, . . . , L max - Tint denotes the integrated autocorre- 
lation time. 



Table |y| displays the parameters of the simulations 
and the integrated autocorrelation times T- mt of the action 
density in multiples of 25 sweeps. It has been determined 
on the L = 8 lattice at (3 < 0.6436 and on L = 12 
at 0.6436 < (3. The integrated autocorrelation of the 
monopole condensate is compatible with these values. 

The f3 values are chosen far enough from the phase 
transition, so that the flips between both phases at the 
phase transition do not occur in the runs. On the largest 
lattices we could simulate this leads to the exclusion of 
runs at j3 = 0.6438. 

Because of the broken Z 2 symmetry, the finite sys- 
tem sometimes flips between the two possible ground 
states within the confinement phase. The flip probabil- 
ity increases with decreasing lattice volume and when 
the phase transition is approached. We observe flips at 
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FIG. 7. Dependence of the monopole condensate on L at 
/3 = 0.64. The curve corresponds to ([L.5|). The scale of the 
y-axis is very fine. 



FIG. 8. Dependence of the monopole condensate on L at 
/3 = 0.6436 (close to the phase transition). The horizontal 
lines indicate the estimated value poo and its error. 



(3 < 0.643 only on the small lattices with L < 8. At 
f3 = 0.6435, we observe them up to the L — 14 and at 
(3 = 0.6436, 0.6437 on all lattice sizes up to L = 18. On 
our largest lattices, the system does not flip more than 
one or two times during the whole run. In these cases we 
cut the corresponding parts of the runs to ensure that the 
unwanted intermediate states do not contribute to the re- 
sulting value of the monopole condensate. However, runs 
with more frequent flips are discarded. 



B. Monopole condensate 

The monopole conde nsat e (0|$+|0) is measured di- 
rectly, and its modulus ( |3.2| ) can be det ermined from the 
correlation function C+(t) of <f> + (3.1). Both methods 
yield compatible results. 



The value of p does not depend on the lattice extent T 
in time direction, but it is weakly dependent on its spacial 
extent L. The L-dependence of p can be described with 
an exponential law 



p(L) = Poo +■ 



-bL 



(5.1) 



with constants a, b and the infinite volume value poa. 
Figure displays p{L) at = 0.64 with the fit ( |5~T| ) used 
for the extrapolation. 

At /3 = 0.6436 and = 0.6437 the data for the con- 
densate could be obtained only on large lattices due to 
the ground state flips on smaller ones. As there is no 
significant L-dependence of p on the largest lattices, we 
take the average of the obt aine d values to represent poo 
instead of using the ansatz (5T). Figure || shows p(L) at 
(3 = 0.6436 and the resulting value of p^. 

Table |v| lists the extrapolated monopole condensate 
values poo at different values of (3. 



p 


poo 


0.6 


0.823(4) 


0.62 


0.730(8) 


0.635 


0.602(5) 


0.638 


0.560(3) 


0.64 


0.509(2) 


0.641 


0.480(6) 


0.642 


0.441(4) 


0.6425 


0.414(5) 


0.643 


0.383(3) 


0.6435 


0.325(10) 


0.6436 


0.300(13) 


0.6437 


0.275(10) 



TABLE V. Extrapolation of the monopole condensate p 
to the infinite spacial volume. 



C. Scaling of the monopole condensate 

The /3-dependence of the monopole condensate poo is 
described with an astonishing precision by a simple power 
law 



Poo (0) = aM™- P) (5-2) 

with the value ( |l.6| ) of the magnetic exponent and 

(3 c c ont = 0.6438(1). (5.3) 

Figure |^ shows the scaling of the monopole condensate 
with the power law (5.2). The logarithmic plot, Fig. nu, 
is obtained by plotting \og(p((3) / a p ) versus log(/3 — (3^ ) 



with a and (3^° nl taken as the fit results from (5.2) 
We note that the value fO) of f3 c c ont and the value fDf) 



of /3^ oul are consistent within two error bars. The value 
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There have been attempts to determine the proper- 
ties of such a state [Q. However, since we use different 
boundary conditions than the authors of |T^ |, and the 
remaining symmetry that is dynamically broken is Z 2 in 
our case instead of the magnetic U(l), our results are 
difficult to compare with those presented in 

We examine the dependence of the energy E± and 
amplitude squares a± = | (1| $± |0) | 2 on the finite spacial 
lattice size L. The corresponding masses m± (L) obtained 
from E± by means of the dispersion relations (2.21) and 
( p. 22 ) show no clear L dependence that would allow an 
extrapolation to the infinite volume. Also the masses m + 
and m_ do not approach each other as L is increased. 



FIG. 9. Scaling behavior of the monopole condensate ex- 
trapolated to L = oo with /3. The curve corresponds to the 
power law (5.2). 




log(/?-/? c ) 

FIG. 10. Scaling of the monopole condensate extrapo- 
lated to L = oo with p in a logarithmic diagram. The line 
corresponds to the power law (5.2). 



of /3™ nt is much more precise, because the determination 
of the condensate in the confinement phase is much more 
precise than that of the monopole mass in the Coulomb 
phase. 



D. Excited states of the condensate 



From the correlation function (3.1) in the confinement 
phase, one can determine in addition to the monopole 
condensate the energy and amplitude of a first excited 
state. If these values gave reasonable results in the in- 
finite volume, they would correspond to a magnetically 
charged particle- like excitation of the monopole conden- 
sate. 



m+(L) 






L 






P 


10 


12 


14 


16 


18 


0.6 


1.90(3) 


1.95(2) 








0.62 


1.43(2) 


1.49(2) 


1.53(2) 






0.635 


0.93(1) 


0.94(1) 


0.98(2) 






0.638 


0.77(1) 


0.78(1) 








0.64 


0.63(2) 


0.635(8) 


0.67(1) 






0.641 


0.554(8) 


0.565(7) 


0.582(7) 


0.61(1) 




0.642 


0.461(7) 


0.479(8) 




0.50(1) 




0.6425 


0.396(8) 


0.439(5) 


0.440(7) 


0.44(1) 


0.446(3) 


0.643 


0.31(1) 


0.29(2) 


0.34(1) 






0.6435 


0.21(2) 


0.20(3) 


0.18(2) 


0.25(3) 


0.26(2) 


0.6436 






0.16(2) 


0.20(2) 


0.19(2) 


0.6437 






0.16(2) 


0.14(2) 


0.16(2) 



TABLE VI. Mass m+ 
confinement phase. 



of the even excited state in the 



m-(L) 






L 






P 


10 


12 


14 


16 


18 


0.6 


2.11(2) 


2.15(2) 








0.62 


1.70(1) 


1.68(2) 


1.69(2) 






0.635 


1.23(1) 


12.3(1) 


1.29(1) 






0.638 


1.08(1) 


1.09(1) 








0.64 


0.99(1) 


0.98(1) 


0.96(1) 






0.641 


0.906(7) 


0.885(6) 


0.890(4) 


0.864(5) 




0.642 


0.810(7) 


0.801(9) 




0.781(9) 




0.6425 


0.764(6) 


0.762(4) 


0.731(8) 


0.716(6) 


0.725(3) 


0.643 


0.719(9) 


0.66(1) 


0.66(1) 






0.6435 


0.58(1) 


0.60(2) 


0.51(2) 


0.52(1) 


0.48(3) 


0.6436 






0.50(1) 


0.50(2) 


0.48(2) 


0.6437 






0.46(2) 


0.45(2) 


0.42(4) 



TABLE VII. Mass m_ 
confinement phase. 



of the odd excited state in the 



Tables VI and VII show the masses m + and m_, re- 



10 




B 9 10 11 12 13 14 15 16 

L 

FIG. 11. Amplitudes a±(L) of the <E>+ (solid) and $- 
(dashed) excitations depending on L with the curves (5.4). 



spectively, for various (3 and the spacial lattice sizes L. 

The amplitudes a±(L) decrease as the spacial lattice 
size is increased. We find that a power law 



a ± {L) = a^ + ^r 



(5.4) 



gives values of the exponent r in the range 1 ... 3 and 
values consistent with zero. Fig. [ll] shows as an 
example the amplitudes a±(L) at j3 = 0.64 with fits by 
the functions (5.4). 

These results indicate, that the observed excitations 
are not present in the infinite volume limit and the corre- 
sponding quantities m± should not be interpreted phys- 
ically as particle masses. For this reason the question of 
an appropriate renormalization of the monopole conden- 
sate remains unanswered. 



VI. CONCLUSION AND DISCUSSION 

In the pure U(l) gauge theory with the Villain 
action we have investigated the scaling behaviour of 
the monopole mass in the Coulomb phase and of the 
monopole condensate in the confinement phase. Both 
observables indicate a critical behavior in the vicinity 
of the phase transition between these phases, with the 



values of the corresponding exponents (1.2) and (1.6). 
Assuming that a continuum theory can be constructed 
at this phase transition, these results indicate that the 
monopoles appear also in such a theory. In particular, 
the monopole mass in the Coulomb phase can vanish 
or stay finite nonzero in physical units. The Gaussian 
value (1.2) of the exponent suggests (though not implies) 
that the corresponding continuum theory may be trivial. 
This property would then hold also for the scalar QED 



in the Coulomb phase, as this theory is obtained by du- 
ality transformations from the pure compact U(l) theory 

00- 

As for the question of the existence of a continuum 
limit, in this work we do not contribute to the resolution 
of the controversy whether the phase transition is of the 
second or of the weak first order [jjj. This was also not 
our aim in this paper, because for this purpose differ- 
ent methods would be more appropriate. As the values 
73f) and (5.3) of (3 C obtained from the Coulomb and the 



confinement phases are consistent within the error bars 
times 1.5, our results are consistent with the second or- 
der. However, a small difference \/3 c c ont - /^? oul | « 0.001 
is not excluded by our data, allowing a weak first order 
transition for the Villain action. 

What we want to point out is that, in spite of such 
a possibility, which actually never can be excluded with 
full certainty, the pure compact QED on the lattice re- 
mains to be a candidate for the construction of an inter- 
esting quantum field theory in continuum. The scaling 
behaviour of the monopole mass and in particular of the 
monopole condensate is of a high quality, described by 
a single exponent in the whole accessible region of val- 
ues of the observables. The same holds also for various 
other observables 0-|j|. This suggests the existence of 
a continuous phase transition somewhere in the parame- 
ter space of possible lattice versions of the pure compact 
QED. This point of view might be relevant also for the 
compact lattice QED with fermions p2). Therefore the 
compact QED merits further investigation with larger re- 
sources and improved methods. 
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